Dynamic exponent in Extremal models of pinning. 
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The depinning transition of a front moving in a time-independent random potential is studied. 
The temporal development of the overall roughness w(L, t) of an initially flat front, w(t) oc t 13 , is the 
classical means to have access to the dynamic exponent. However, in the case of front propagation in 
quenched disorder via extremal dynamics, we show that the initial increase in front roughness implies 
an extra dependence over the system size which comes from the fact that the activity is essentially 
localized in a narrow region of space. We propose an analytic expression for the j3 exponent and 
confirm this for different models (crack front propagation, Edwards-Wilkinson model in a quenched 



The propagation of a front in a noisy environment has been the subject of active research in the past. In particular, 
after Family and Vicsek proposed a scaling form for the evolution of the roughness of fronts and surfaces from Langevin 
equations, numerous analytical and numerical works have verified these laws in a wide variety of models. Various 
reviews cover this rich field . 

More recently, the quenched (i.e. time-independent) nature of the noise was recognised as playing a significant role 
in front propagation. Unfortunately, in spite of a few key works, analytic modelling of such a depinning transition 
is rather scarce: Dynamic Renormalization Group studies have been proposed for the Edwards- Wilkinson model 
with quenched disorder [Q-^) however in 1+1 dimensions some of the predicted exponents are quite far from their 
estimate obtained in numerical simulations. It has been proposed that the Kardar-Parisi-Zhang model in a quenched 
environment can be described by a directed percolation related model, mostly on the basis of numerical agreement 
between measured exponents 0|| . 

Such models are relevant for a number of pinning phenomena such as crack propagation [p|-|l3|] , wetting phenomena 
]r^ . [T5[ , vortex pinning in type-II superconductors jl6| and solid friction |l7],|l8| • 

Growth models in 1+1 dimensions are often studied through the evolution in time 9, (defined as the number of 
growth steps t divided by the system size L), of the roughness, w{6), of an initially flat front u>(0) = 0. The roughness 
— standard deviation of the front — in a system of size L obeys the scaling form 

where 

, s (est for x <C 1 /rA 
^ = U< for*»l (2) 

and £ = pzi. The exponent z\ is referred to as the dynamic exponent since it relates space and time, whereas Q 
describes the roughness of the front (Hurst or roughness exponent), i.e. the scaling of the pair correlation function 
for equal-time positions along front. In most cases of annealed noise, there is no need Jh| to introduce any other 
exponents, since at late stages, when the overall roughness has reached the saturation value, the full two point 
correlation function of the front at different locations and times reveals a similar scaling between space and time: 
Aa; cx A6* 1 / 21 . 

The aim of this letter is to show that for a class of quenched disorder depinning models, the early time development 
of the roughness does not obey such a law. Rather, it implies an extra L dependence in the expression of w, from 
which £ = (3zi is violated. 

We focus more specifically on a class of models introduced by Tanguy et al pl[ obeying extremal dynamics. The 
front is defined by its position z = h(x,t). An external driving F allows to exert a pressure on the front which is 
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however biased by the position of the front, so that at each site x, the force is f(x, t) = F(G*h) where the * denotes 
a convolution product, and G is a function which is specific to the physical problem studied. The local part of G is 
adjusted so that G(0) = — X^/o G(x). In order to normalize this kernel, we choose as a convention G(0) = — 1. The 
key feature is that G may have power-law tails which determine the universality class of the problem. For a planar 
crack, Gao and Rice Jl3| have shown that (to first order in h), G decays as r~ 2 . The same holds for the motion of 
the triple line of a liquid-glass interface intersecting a solid surface, in a wetting problem with a semi-infinite liquid 
surface in weak gravity. If the liquid/gas system is confined between two parallel plates the G function decays more 
abruptly, as r -3 . In the mean-field limit, the G function does not depend on distance. For the sake of convenience, 
this model has thus been generalized to any power-law form for G: 

G(r) oc r- a (3) 

In our simulations, periodic boundary conditions are implemented and thus we adjust G to match such a periodicity: 
namely G(r) oc sin(7rr/L) _Q . Such an expression is exact for a — 2, and is a convenient means to implement the 
periodic b.c. in other cases. 

The environment is represented by quenched heterogeneities, r)(x, z), which can block the front whenever f(x,t) < 
r){x,h{x,t)). The r\ values are considered to be uncorrelated, positive numbers, randomly picked from a uniform 
distribution over the interval [0; 1]. Thus for F = 0, the interface does not move. At each time step, the force is 
increased slowly from up to the level where one site x* (the active site) can depin. This site jumps to the next 
obstacle r\ at a random distance along the z direction. This distance is chosen again from a uniform distribution 
over the interval [0;g?]. The external loading is immediately brought back to zero so that no other sites can move 
simultaneously. The same step is repeated indefinitely, d is a free parameter. It has been checked that d plays no role 
in the statistical properties of the model in the steady state. 

The steady state properties of this model have been studied in detail numerically pp|-|2^]. It has been shown in 
particular that for a < 1, the model is in the mean field regime, as can be easily inferred from the a = case. 
No spatial structure appears, and the front is an uncorrelated white noise. For a > 3, the long-range kernel is 
dominated by the short wavelength cut-off, and is thus equivalent to the case where G is the second derivative of a 
Dirac distribution, i.e. the local force is proportional to the external loading and the local curvature of the front. 
This is the Edwards- Wilkinson or "Laplacian" case with quenched disorder. The front has a "super-rough" structure 
with a roughness exponent (a 1,2 ^3|. The steady state properties of this dynamics is fairly rich and a number of 
scaling properties can be observed for the front structure, the time evolution of the activity, and of the driving force 
necessary to depin the interface. In the intermediate range 1 < a < 3, similar properties are observed with scaling 
exponents which vary continuously with a. 

One especially interesting property can be studied in order to characterize the spreading of activity in space and 
time: Let x*(to) be the active site at time to, and x*(t + At) at time to + At. Note that here we use as a time the 
total number of moves t rather than the number of moves per site 9 = t/L. We study the statistical distribution 
p(Ax, At) of Ax = \x*(t + At) - x*(t Q )\ for a fixed At. It obeys the scaling form 



,(A*, At) = At- 1 /*- 



(4) 



where 

v ' l x~ b for x > 1 

where b is equal to a. Such a scaling form was first introduced by Furuberg et al studying invasion percolation 
Thus, in the steady state, we see that the activity spreads typically over distances Ax oc At 1 / 22 . In this sense, zi is 
actually the "dynamic exponent" , and it indeed governs all correspondences between time and space in the steady 
state. Note however that a different convention is now used for z, because of the definition of time. In particular the 
activity has spread over the entire system for a time equal to t oc L Z2 , hence 9 = t/L oc L Z2_1 . Thus actually Z2 should 
be compared to Z\ + 1. In models with extremal dynamics, can easily be related to the roughness exponent of the 
front [ p5|j2l| . After a time At, the activity remains localized in a region of extend Ax. Over this region, the front 
moves by a typical distance of order Az oc Ax^ . Thus the number of time steps required to travel by this amount 
scales as AxAz oc Ax 1+I ^ , hence 

Z2 = 1 + C (6) 
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This simple argument has been checked to be obeyed precisely in numerical simulations for 1 < a < 3. Let us however 
note that it breaks down in the Laplacian case a — > oo where zi ~ 2.0 and £ ~ 1.25. However, the situation when 
C > 1 is known to display some pathological behaviors. 
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FIG. 1. Log-log plot of the overall roughness of the interface as a function of time for two values of the a parameter a = 2 
(symbol o) and a — 3 (symbol •). The best power- law fits are shown as plain lines. The values of the slopes are reported in 
Table |l| The size of the system is L=1024. 

We now consider the early stages of the front growth. We assign an initial flat configuration for the front z = 
and study the evolution of the front roughness, in order to characterize z\ through the scaling relation Eq. |l]. Figure 
1 shows a log-log plot of the overall roughness for the crack model, a = 2, and for a = 3. In both cases, we indeed 
measure a power-law, w oct' 3 , with (3 w 0.61 and 0.70 respectively. Table || gives the different values of C and Z2 from 
Ref. [^l| , and (3 from the numerical simulations of the present study. We note that if we blindly apply the relation 
z i = C/A we nn d a severe discrepancy with Z2 — 1 which cannot be attributed to numerical uncertainties. For instance, 
if a = 2, £/ ' (3 = 0.58, and Z2 — 1 = 0.35. Let us note that this observation invalidates the numerical determination of 
the z exponent published for example in Refs. (JJ^] (moreover, in these articles the analytical expression proposed for 
z , using RG analysis, refers to another driving mode that could be interesting to compare with extremal models). 
This is the main message of this paper: The standard relation z = breaks down for quenched disorder growth 
models with extremal dynamics. 
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FIG. 2. Shape of the first stages t = 1000, 6000 and 11000 of the interface for a = 2 and L = 1024. We note that the activity 
is not equally distributed along the interface, but rather only a localized part grows and reaches a steady state conformation 
before growing and invading the fiat region. This observation expressed in quantitative terms gives our prediction of j3. 



We can go one step beyond the numerical study reported so far and estimate the value of z\ from the steady state 
roughness exponent. To this end, it is informative to look at the shape of the front in the initial regime. Figure 
2 shows such an example for a = 2 and L = 1024. We observe that the front remains pinned along its initial fiat 
geometry for a large time, and thus the interface moves only in a confined region of space. This region progressively 
grows, and invades the flat part, but still continues its propagation along the z direction. This is in marked contrast 
with annealed noise type growth where the activity is delocalized along the entire front. Moreover, as can be noted 
from Figure 1 , the cross-over to the saturation regime is quite steep. This suggests that the depinned parts of the 
front have already reached their steady state roughness. 

Let us now translate this argument into quantitative terms. We assume that the interface is depinned along one 
single interval of length i. Along this interval, the roughness assumes its asymptotic value, i.e. the typical height will 
be of order h « tS. Therefore the time needed to reach this position is 



t oc £ 1+ < 



(7) 



Generally, as one follows the time development of the roughness, one introduces as a time the number of moves per 
site, 9 = t/L. Ignoring subdominant terms coming from the average height of the interface, we compute the overall 
roughness, w, as 



w 2 oc 



£2C+1 



(8) 



where the multiplicative £/L term comes from the weight of the depinned part of the interface as compared to the 
rest. Therefore the [3 exponent is readily estimated from the elimination of t in the two above equations, and thus 



with 



w oc O^L^ 1 ' 2 
C + l/2 



ft 



C + i 



(9) 
(10) 
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We note however that there is an extra term depending on L in the rate of growth of the roughness. Thus one cannot 
use the standard relation z\ = Q/ fl. Instead, we write 



w 



L 



(11) 



The dynamic exponent z\ is then found by imposing that in the steady state w is i-independent, and w oc . This 
gives the power-law behavior of ip for small arguments, f(x) oc x^ +1 / 2 ~P , and thus 



(3 



C + 1/2 - 



Zl 







or using our above expression for (3 



C + 1/2-/3 , 



(12) 



(13) 



Thus the total number of moves necessary to reach the steady state scales as t = 9L oc L 1+< % and thus we recover our 
previous expression for Z2 = 1 + z± where the difference of one simply comes from the definition of time. 



a 


from Ref. plj 


Z2 

from Ref. fl| 


a 

(measured) 


from Eq. © 


1.5 


0.05 


1.05 


0.49 


0.52 


2.0 


0.35 


1.35 


0.61 


0.63 


2.5 


0.65 


1.65 


0.65 


0.70 


3.0 


1.0 


2.0 


0.70 


0.75 


oo 


1.25 


2.0 


0.8* 


0.78 



TABLE I. Values of the exponents £ and Z2 taken from Ref. [21], measured value of j3 from the present work, and predicted 
value from Eq. (10). * refers to Ref. [23]. 
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In Table Q, we have reported the value of the (3 exponent estimated from numerical simulation data, and the one 
obtained from Eq. [l0| We note an excellent numerical agreement for all values of a studied. 

More generally if other moments of the interface height are computed, a different scaling is expected. Indeed for 
the moment of order m, a similar computation gives ((h- {h)) m ) l / m oct?™ with 

A» = ^ (14) 

Hence, in contrast with annealed models, different moments give rise to different estimates of the (3 exponent. This 
again underlines the fact that care has to be taken with the interpretation of the latter exponent. 

We mentioned that the part of the interface which has moved had already a conformation representative of the 
steady state regime. In order to check this, we can perform the following test. After the interface has reached the 
steady state, we choose an arbitrary time to and record the position of the interface ho(x) — h(x,to). Then as the 
simulation continues, we study the incremental motion of the front, Ah(x, t) — h(x, t) — ho(x). The time evolution of 
the roughness of Ah follows indeed the same law as the early stage i /3 . This shows that the initial flat configuration 
behaves as any late stage configuration. 

The final question to answer is why the activity is localized during the early stages of front propagation. As a site 
moves, a large part of the forces it carries is transferred to the nearest neighbours, while a smaller part is transferred to 
the second neighbours, and so on. The amount of load transfer depends on the distance to the active site as dictated 
by the kernel G of the model, and hence it depends on a. For a — > oo, the load transfer is local, and only the nearest 
neighbors are influenced. The amplitude of the change in force is by the definition of the model dependent on the 
distance the active site advances, and thus it depends on the parameter d. For large d, one may easily understand that 
the activity has a tendency to move uniformly in the transverse direction. The first active site will jump by a distance 
proportional to d, and hence its closest neighbours will be pushed forward by such an amplitude that the amplitude 
of the threshold strength may be insufficient to keep them pinned. In contrast, for small d, the threshold distribution 
may win over the force modification induced by the roughening of the front. As the amplitude of d is reduced, indeed, 
the activity map shows that the very initial stage is spread over the entire interface. In the very early stages, one 
measures a roughness exponent (3 « 0.5 as can be expected from the trivial observation that only a number t of sites 
move by a single step proportional to d, thus w oc dt 1 / 2 L~ 1 / 2 . After this initial transient, a higher slope takes over. 
Figure 3 || indicates the evolution of w(t) for three values of d, 1.0, 0.1 and 0.01. Fitting the time region where on 
average there has been more than one move per site, t > L, gives consistent estimates of (3 as mentioned in Table |. 
We thus conclude that the above description of the early stage growth is the generic case, which is encountered for 
any specific choice of the parameter d. 
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FIG. 3. Log-log plot of the overall roughness versus time for a = 2 and L = 1024. Three values of d are used as indicated 
in the caption. The dotted line shows a slope of 0.5 which accounts for the very early stages of growth when d is small. 

Let us conclude by briefly summarizing our results: The time development of roughness in quenched disorder 
depinning models implies a size dependence which has not been noted before. This implies a violation of the scaling 
relation £ = z(3. Instead, we show that the activity is localized even in the early stages of growth and thus this implies 
a power-law increase of the roughness in time with an exponent /3 given in Eq. ([l^) . This has been confirmed through 
numerical simulations. 
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